3.2.5 \(\int (a+b x+c x^2)^{3/2} (d+e x+f x^2) \, dx\) [105]

3.2.5.1 Optimal result
3.2.5.2 Mathematica [A] (verified)
3.2.5.3 Rubi [A] (verified)
3.2.5.4 Maple [A] (verified)
3.2.5.5 Fricas [A] (verification not implemented)
3.2.5.6 Sympy [B] (verification not implemented)
3.2.5.7 Maxima [F(-2)]
3.2.5.8 Giac [A] (verification not implemented)
3.2.5.9 Mupad [F(-1)]

3.2.5.1 Optimal result

Integrand size = 25, antiderivative size = 236 \[ \int \left (a+b x+c x^2\right )^{3/2} \left (d+e x+f x^2\right ) \, dx=-\frac {\left (b^2-4 a c\right ) \left (24 c^2 d+7 b^2 f-4 c (3 b e+a f)\right ) (b+2 c x) \sqrt {a+b x+c x^2}}{512 c^4}+\frac {\left (24 c^2 d-12 b c e+7 b^2 f-4 a c f\right ) (b+2 c x) \left (a+b x+c x^2\right )^{3/2}}{192 c^3}+\frac {(12 c e-7 b f) \left (a+b x+c x^2\right )^{5/2}}{60 c^2}+\frac {f x \left (a+b x+c x^2\right )^{5/2}}{6 c}+\frac {\left (b^2-4 a c\right )^2 \left (24 c^2 d+7 b^2 f-4 c (3 b e+a f)\right ) \text {arctanh}\left (\frac {b+2 c x}{2 \sqrt {c} \sqrt {a+b x+c x^2}}\right )}{1024 c^{9/2}} \]

output
1/192*(-4*a*c*f+7*b^2*f-12*b*c*e+24*c^2*d)*(2*c*x+b)*(c*x^2+b*x+a)^(3/2)/c 
^3+1/60*(-7*b*f+12*c*e)*(c*x^2+b*x+a)^(5/2)/c^2+1/6*f*x*(c*x^2+b*x+a)^(5/2 
)/c+1/1024*(-4*a*c+b^2)^2*(24*c^2*d+7*b^2*f-4*c*(a*f+3*b*e))*arctanh(1/2*( 
2*c*x+b)/c^(1/2)/(c*x^2+b*x+a)^(1/2))/c^(9/2)-1/512*(-4*a*c+b^2)*(24*c^2*d 
+7*b^2*f-4*c*(a*f+3*b*e))*(2*c*x+b)*(c*x^2+b*x+a)^(1/2)/c^4
 
3.2.5.2 Mathematica [A] (verified)

Time = 2.72 (sec) , antiderivative size = 293, normalized size of antiderivative = 1.24 \[ \int \left (a+b x+c x^2\right )^{3/2} \left (d+e x+f x^2\right ) \, dx=\frac {\sqrt {c} \sqrt {a+x (b+c x)} \left (-105 b^5 f+10 b^4 c (18 e+7 f x)-8 b^3 c (45 c d-95 a f+c x (15 e+7 f x))+48 b^2 c^2 (-a (25 e+9 f x)+c x (5 d+x (2 e+f x)))+16 b c^2 \left (-81 a^2 f+6 a c (25 d+x (7 e+3 f x))+4 c^2 x^2 (45 d+x (33 e+26 f x))\right )+32 c^3 \left (3 a^2 (16 e+5 f x)+4 c^2 x^3 (15 d+2 x (6 e+5 f x))+2 a c x (75 d+x (48 e+35 f x))\right )\right )+15 \left (b^2-4 a c\right )^2 \left (24 c^2 d+7 b^2 f-4 c (3 b e+a f)\right ) \text {arctanh}\left (\frac {\sqrt {c} x}{-\sqrt {a}+\sqrt {a+x (b+c x)}}\right )}{7680 c^{9/2}} \]

input
Integrate[(a + b*x + c*x^2)^(3/2)*(d + e*x + f*x^2),x]
 
output
(Sqrt[c]*Sqrt[a + x*(b + c*x)]*(-105*b^5*f + 10*b^4*c*(18*e + 7*f*x) - 8*b 
^3*c*(45*c*d - 95*a*f + c*x*(15*e + 7*f*x)) + 48*b^2*c^2*(-(a*(25*e + 9*f* 
x)) + c*x*(5*d + x*(2*e + f*x))) + 16*b*c^2*(-81*a^2*f + 6*a*c*(25*d + x*( 
7*e + 3*f*x)) + 4*c^2*x^2*(45*d + x*(33*e + 26*f*x))) + 32*c^3*(3*a^2*(16* 
e + 5*f*x) + 4*c^2*x^3*(15*d + 2*x*(6*e + 5*f*x)) + 2*a*c*x*(75*d + x*(48* 
e + 35*f*x)))) + 15*(b^2 - 4*a*c)^2*(24*c^2*d + 7*b^2*f - 4*c*(3*b*e + a*f 
))*ArcTanh[(Sqrt[c]*x)/(-Sqrt[a] + Sqrt[a + x*(b + c*x)])])/(7680*c^(9/2))
 
3.2.5.3 Rubi [A] (verified)

Time = 0.36 (sec) , antiderivative size = 210, normalized size of antiderivative = 0.89, number of steps used = 8, number of rules used = 7, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.280, Rules used = {2192, 27, 1160, 1087, 1087, 1092, 219}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int \left (a+b x+c x^2\right )^{3/2} \left (d+e x+f x^2\right ) \, dx\)

\(\Big \downarrow \) 2192

\(\displaystyle \frac {\int \frac {1}{2} (12 c d-2 a f+(12 c e-7 b f) x) \left (c x^2+b x+a\right )^{3/2}dx}{6 c}+\frac {f x \left (a+b x+c x^2\right )^{5/2}}{6 c}\)

\(\Big \downarrow \) 27

\(\displaystyle \frac {\int (2 (6 c d-a f)+(12 c e-7 b f) x) \left (c x^2+b x+a\right )^{3/2}dx}{12 c}+\frac {f x \left (a+b x+c x^2\right )^{5/2}}{6 c}\)

\(\Big \downarrow \) 1160

\(\displaystyle \frac {\frac {\left (-4 a c f+7 b^2 f-12 b c e+24 c^2 d\right ) \int \left (c x^2+b x+a\right )^{3/2}dx}{2 c}+\frac {\left (a+b x+c x^2\right )^{5/2} (12 c e-7 b f)}{5 c}}{12 c}+\frac {f x \left (a+b x+c x^2\right )^{5/2}}{6 c}\)

\(\Big \downarrow \) 1087

\(\displaystyle \frac {\frac {\left (-4 a c f+7 b^2 f-12 b c e+24 c^2 d\right ) \left (\frac {(b+2 c x) \left (a+b x+c x^2\right )^{3/2}}{8 c}-\frac {3 \left (b^2-4 a c\right ) \int \sqrt {c x^2+b x+a}dx}{16 c}\right )}{2 c}+\frac {\left (a+b x+c x^2\right )^{5/2} (12 c e-7 b f)}{5 c}}{12 c}+\frac {f x \left (a+b x+c x^2\right )^{5/2}}{6 c}\)

\(\Big \downarrow \) 1087

\(\displaystyle \frac {\frac {\left (-4 a c f+7 b^2 f-12 b c e+24 c^2 d\right ) \left (\frac {(b+2 c x) \left (a+b x+c x^2\right )^{3/2}}{8 c}-\frac {3 \left (b^2-4 a c\right ) \left (\frac {(b+2 c x) \sqrt {a+b x+c x^2}}{4 c}-\frac {\left (b^2-4 a c\right ) \int \frac {1}{\sqrt {c x^2+b x+a}}dx}{8 c}\right )}{16 c}\right )}{2 c}+\frac {\left (a+b x+c x^2\right )^{5/2} (12 c e-7 b f)}{5 c}}{12 c}+\frac {f x \left (a+b x+c x^2\right )^{5/2}}{6 c}\)

\(\Big \downarrow \) 1092

\(\displaystyle \frac {\frac {\left (-4 a c f+7 b^2 f-12 b c e+24 c^2 d\right ) \left (\frac {(b+2 c x) \left (a+b x+c x^2\right )^{3/2}}{8 c}-\frac {3 \left (b^2-4 a c\right ) \left (\frac {(b+2 c x) \sqrt {a+b x+c x^2}}{4 c}-\frac {\left (b^2-4 a c\right ) \int \frac {1}{4 c-\frac {(b+2 c x)^2}{c x^2+b x+a}}d\frac {b+2 c x}{\sqrt {c x^2+b x+a}}}{4 c}\right )}{16 c}\right )}{2 c}+\frac {\left (a+b x+c x^2\right )^{5/2} (12 c e-7 b f)}{5 c}}{12 c}+\frac {f x \left (a+b x+c x^2\right )^{5/2}}{6 c}\)

\(\Big \downarrow \) 219

\(\displaystyle \frac {\frac {\left (\frac {(b+2 c x) \left (a+b x+c x^2\right )^{3/2}}{8 c}-\frac {3 \left (b^2-4 a c\right ) \left (\frac {(b+2 c x) \sqrt {a+b x+c x^2}}{4 c}-\frac {\left (b^2-4 a c\right ) \text {arctanh}\left (\frac {b+2 c x}{2 \sqrt {c} \sqrt {a+b x+c x^2}}\right )}{8 c^{3/2}}\right )}{16 c}\right ) \left (-4 a c f+7 b^2 f-12 b c e+24 c^2 d\right )}{2 c}+\frac {\left (a+b x+c x^2\right )^{5/2} (12 c e-7 b f)}{5 c}}{12 c}+\frac {f x \left (a+b x+c x^2\right )^{5/2}}{6 c}\)

input
Int[(a + b*x + c*x^2)^(3/2)*(d + e*x + f*x^2),x]
 
output
(f*x*(a + b*x + c*x^2)^(5/2))/(6*c) + (((12*c*e - 7*b*f)*(a + b*x + c*x^2) 
^(5/2))/(5*c) + ((24*c^2*d - 12*b*c*e + 7*b^2*f - 4*a*c*f)*(((b + 2*c*x)*( 
a + b*x + c*x^2)^(3/2))/(8*c) - (3*(b^2 - 4*a*c)*(((b + 2*c*x)*Sqrt[a + b* 
x + c*x^2])/(4*c) - ((b^2 - 4*a*c)*ArcTanh[(b + 2*c*x)/(2*Sqrt[c]*Sqrt[a + 
 b*x + c*x^2])])/(8*c^(3/2))))/(16*c)))/(2*c))/(12*c)
 

3.2.5.3.1 Defintions of rubi rules used

rule 27
Int[(a_)*(Fx_), x_Symbol] :> Simp[a   Int[Fx, x], x] /; FreeQ[a, x] &&  !Ma 
tchQ[Fx, (b_)*(Gx_) /; FreeQ[b, x]]
 

rule 219
Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(1/(Rt[a, 2]*Rt[-b, 2]))* 
ArcTanh[Rt[-b, 2]*(x/Rt[a, 2])], x] /; FreeQ[{a, b}, x] && NegQ[a/b] && (Gt 
Q[a, 0] || LtQ[b, 0])
 

rule 1087
Int[((a_.) + (b_.)*(x_) + (c_.)*(x_)^2)^(p_), x_Symbol] :> Simp[(b + 2*c*x) 
*((a + b*x + c*x^2)^p/(2*c*(2*p + 1))), x] - Simp[p*((b^2 - 4*a*c)/(2*c*(2* 
p + 1)))   Int[(a + b*x + c*x^2)^(p - 1), x], x] /; FreeQ[{a, b, c}, x] && 
GtQ[p, 0] && (IntegerQ[4*p] || IntegerQ[3*p])
 

rule 1092
Int[1/Sqrt[(a_) + (b_.)*(x_) + (c_.)*(x_)^2], x_Symbol] :> Simp[2   Subst[I 
nt[1/(4*c - x^2), x], x, (b + 2*c*x)/Sqrt[a + b*x + c*x^2]], x] /; FreeQ[{a 
, b, c}, x]
 

rule 1160
Int[((d_.) + (e_.)*(x_))*((a_.) + (b_.)*(x_) + (c_.)*(x_)^2)^(p_), x_Symbol 
] :> Simp[e*((a + b*x + c*x^2)^(p + 1)/(2*c*(p + 1))), x] + Simp[(2*c*d - b 
*e)/(2*c)   Int[(a + b*x + c*x^2)^p, x], x] /; FreeQ[{a, b, c, d, e, p}, x] 
 && NeQ[p, -1]
 

rule 2192
Int[(Pq_)*((a_.) + (b_.)*(x_) + (c_.)*(x_)^2)^(p_), x_Symbol] :> With[{q = 
Expon[Pq, x], e = Coeff[Pq, x, Expon[Pq, x]]}, Simp[e*x^(q - 1)*((a + b*x + 
 c*x^2)^(p + 1)/(c*(q + 2*p + 1))), x] + Simp[1/(c*(q + 2*p + 1))   Int[(a 
+ b*x + c*x^2)^p*ExpandToSum[c*(q + 2*p + 1)*Pq - a*e*(q - 1)*x^(q - 2) - b 
*e*(q + p)*x^(q - 1) - c*e*(q + 2*p + 1)*x^q, x], x], x]] /; FreeQ[{a, b, c 
, p}, x] && PolyQ[Pq, x] && NeQ[b^2 - 4*a*c, 0] &&  !LeQ[p, -1]
 
3.2.5.4 Maple [A] (verified)

Time = 0.66 (sec) , antiderivative size = 398, normalized size of antiderivative = 1.69

method result size
risch \(-\frac {\left (-1280 c^{5} f \,x^{5}-1664 b \,c^{4} f \,x^{4}-1536 c^{5} e \,x^{4}-2240 a \,c^{4} f \,x^{3}-48 b^{2} c^{3} f \,x^{3}-2112 b \,c^{4} e \,x^{3}-1920 c^{5} d \,x^{3}-288 a b \,c^{3} f \,x^{2}-3072 a \,c^{4} e \,x^{2}+56 b^{3} c^{2} f \,x^{2}-96 b^{2} c^{3} e \,x^{2}-2880 b \,c^{4} d \,x^{2}-480 a^{2} c^{3} f x +432 a \,b^{2} c^{2} f x -672 a b \,c^{3} e x -4800 a \,c^{4} d x -70 b^{4} c f x +120 b^{3} c^{2} e x -240 b^{2} c^{3} d x +1296 a^{2} b \,c^{2} f -1536 a^{2} c^{3} e -760 a \,b^{3} c f +1200 a \,b^{2} c^{2} e -2400 a b \,c^{3} d +105 b^{5} f -180 b^{4} c e +360 b^{3} c^{2} d \right ) \sqrt {c \,x^{2}+b x +a}}{7680 c^{4}}-\frac {\left (64 a^{3} c^{3} f -144 a^{2} b^{2} c^{2} f +192 a^{2} b \,c^{3} e -384 a^{2} c^{4} d +60 a \,b^{4} c f -96 a \,b^{3} c^{2} e +192 a \,b^{2} c^{3} d -7 b^{6} f +12 b^{5} c e -24 b^{4} c^{2} d \right ) \ln \left (\frac {\frac {b}{2}+c x}{\sqrt {c}}+\sqrt {c \,x^{2}+b x +a}\right )}{1024 c^{\frac {9}{2}}}\) \(398\)
default \(d \left (\frac {\left (2 c x +b \right ) \left (c \,x^{2}+b x +a \right )^{\frac {3}{2}}}{8 c}+\frac {3 \left (4 a c -b^{2}\right ) \left (\frac {\left (2 c x +b \right ) \sqrt {c \,x^{2}+b x +a}}{4 c}+\frac {\left (4 a c -b^{2}\right ) \ln \left (\frac {\frac {b}{2}+c x}{\sqrt {c}}+\sqrt {c \,x^{2}+b x +a}\right )}{8 c^{\frac {3}{2}}}\right )}{16 c}\right )+f \left (\frac {x \left (c \,x^{2}+b x +a \right )^{\frac {5}{2}}}{6 c}-\frac {7 b \left (\frac {\left (c \,x^{2}+b x +a \right )^{\frac {5}{2}}}{5 c}-\frac {b \left (\frac {\left (2 c x +b \right ) \left (c \,x^{2}+b x +a \right )^{\frac {3}{2}}}{8 c}+\frac {3 \left (4 a c -b^{2}\right ) \left (\frac {\left (2 c x +b \right ) \sqrt {c \,x^{2}+b x +a}}{4 c}+\frac {\left (4 a c -b^{2}\right ) \ln \left (\frac {\frac {b}{2}+c x}{\sqrt {c}}+\sqrt {c \,x^{2}+b x +a}\right )}{8 c^{\frac {3}{2}}}\right )}{16 c}\right )}{2 c}\right )}{12 c}-\frac {a \left (\frac {\left (2 c x +b \right ) \left (c \,x^{2}+b x +a \right )^{\frac {3}{2}}}{8 c}+\frac {3 \left (4 a c -b^{2}\right ) \left (\frac {\left (2 c x +b \right ) \sqrt {c \,x^{2}+b x +a}}{4 c}+\frac {\left (4 a c -b^{2}\right ) \ln \left (\frac {\frac {b}{2}+c x}{\sqrt {c}}+\sqrt {c \,x^{2}+b x +a}\right )}{8 c^{\frac {3}{2}}}\right )}{16 c}\right )}{6 c}\right )+e \left (\frac {\left (c \,x^{2}+b x +a \right )^{\frac {5}{2}}}{5 c}-\frac {b \left (\frac {\left (2 c x +b \right ) \left (c \,x^{2}+b x +a \right )^{\frac {3}{2}}}{8 c}+\frac {3 \left (4 a c -b^{2}\right ) \left (\frac {\left (2 c x +b \right ) \sqrt {c \,x^{2}+b x +a}}{4 c}+\frac {\left (4 a c -b^{2}\right ) \ln \left (\frac {\frac {b}{2}+c x}{\sqrt {c}}+\sqrt {c \,x^{2}+b x +a}\right )}{8 c^{\frac {3}{2}}}\right )}{16 c}\right )}{2 c}\right )\) \(499\)

input
int((c*x^2+b*x+a)^(3/2)*(f*x^2+e*x+d),x,method=_RETURNVERBOSE)
 
output
-1/7680/c^4*(-1280*c^5*f*x^5-1664*b*c^4*f*x^4-1536*c^5*e*x^4-2240*a*c^4*f* 
x^3-48*b^2*c^3*f*x^3-2112*b*c^4*e*x^3-1920*c^5*d*x^3-288*a*b*c^3*f*x^2-307 
2*a*c^4*e*x^2+56*b^3*c^2*f*x^2-96*b^2*c^3*e*x^2-2880*b*c^4*d*x^2-480*a^2*c 
^3*f*x+432*a*b^2*c^2*f*x-672*a*b*c^3*e*x-4800*a*c^4*d*x-70*b^4*c*f*x+120*b 
^3*c^2*e*x-240*b^2*c^3*d*x+1296*a^2*b*c^2*f-1536*a^2*c^3*e-760*a*b^3*c*f+1 
200*a*b^2*c^2*e-2400*a*b*c^3*d+105*b^5*f-180*b^4*c*e+360*b^3*c^2*d)*(c*x^2 
+b*x+a)^(1/2)-1/1024*(64*a^3*c^3*f-144*a^2*b^2*c^2*f+192*a^2*b*c^3*e-384*a 
^2*c^4*d+60*a*b^4*c*f-96*a*b^3*c^2*e+192*a*b^2*c^3*d-7*b^6*f+12*b^5*c*e-24 
*b^4*c^2*d)/c^(9/2)*ln((1/2*b+c*x)/c^(1/2)+(c*x^2+b*x+a)^(1/2))
 
3.2.5.5 Fricas [A] (verification not implemented)

Time = 0.40 (sec) , antiderivative size = 839, normalized size of antiderivative = 3.56 \[ \int \left (a+b x+c x^2\right )^{3/2} \left (d+e x+f x^2\right ) \, dx=\left [-\frac {15 \, {\left (24 \, {\left (b^{4} c^{2} - 8 \, a b^{2} c^{3} + 16 \, a^{2} c^{4}\right )} d - 12 \, {\left (b^{5} c - 8 \, a b^{3} c^{2} + 16 \, a^{2} b c^{3}\right )} e + {\left (7 \, b^{6} - 60 \, a b^{4} c + 144 \, a^{2} b^{2} c^{2} - 64 \, a^{3} c^{3}\right )} f\right )} \sqrt {c} \log \left (-8 \, c^{2} x^{2} - 8 \, b c x - b^{2} + 4 \, \sqrt {c x^{2} + b x + a} {\left (2 \, c x + b\right )} \sqrt {c} - 4 \, a c\right ) - 4 \, {\left (1280 \, c^{6} f x^{5} + 128 \, {\left (12 \, c^{6} e + 13 \, b c^{5} f\right )} x^{4} + 16 \, {\left (120 \, c^{6} d + 132 \, b c^{5} e + {\left (3 \, b^{2} c^{4} + 140 \, a c^{5}\right )} f\right )} x^{3} + 8 \, {\left (360 \, b c^{5} d + 12 \, {\left (b^{2} c^{4} + 32 \, a c^{5}\right )} e - {\left (7 \, b^{3} c^{3} - 36 \, a b c^{4}\right )} f\right )} x^{2} - 120 \, {\left (3 \, b^{3} c^{3} - 20 \, a b c^{4}\right )} d + 12 \, {\left (15 \, b^{4} c^{2} - 100 \, a b^{2} c^{3} + 128 \, a^{2} c^{4}\right )} e - {\left (105 \, b^{5} c - 760 \, a b^{3} c^{2} + 1296 \, a^{2} b c^{3}\right )} f + 2 \, {\left (120 \, {\left (b^{2} c^{4} + 20 \, a c^{5}\right )} d - 12 \, {\left (5 \, b^{3} c^{3} - 28 \, a b c^{4}\right )} e + {\left (35 \, b^{4} c^{2} - 216 \, a b^{2} c^{3} + 240 \, a^{2} c^{4}\right )} f\right )} x\right )} \sqrt {c x^{2} + b x + a}}{30720 \, c^{5}}, -\frac {15 \, {\left (24 \, {\left (b^{4} c^{2} - 8 \, a b^{2} c^{3} + 16 \, a^{2} c^{4}\right )} d - 12 \, {\left (b^{5} c - 8 \, a b^{3} c^{2} + 16 \, a^{2} b c^{3}\right )} e + {\left (7 \, b^{6} - 60 \, a b^{4} c + 144 \, a^{2} b^{2} c^{2} - 64 \, a^{3} c^{3}\right )} f\right )} \sqrt {-c} \arctan \left (\frac {\sqrt {c x^{2} + b x + a} {\left (2 \, c x + b\right )} \sqrt {-c}}{2 \, {\left (c^{2} x^{2} + b c x + a c\right )}}\right ) - 2 \, {\left (1280 \, c^{6} f x^{5} + 128 \, {\left (12 \, c^{6} e + 13 \, b c^{5} f\right )} x^{4} + 16 \, {\left (120 \, c^{6} d + 132 \, b c^{5} e + {\left (3 \, b^{2} c^{4} + 140 \, a c^{5}\right )} f\right )} x^{3} + 8 \, {\left (360 \, b c^{5} d + 12 \, {\left (b^{2} c^{4} + 32 \, a c^{5}\right )} e - {\left (7 \, b^{3} c^{3} - 36 \, a b c^{4}\right )} f\right )} x^{2} - 120 \, {\left (3 \, b^{3} c^{3} - 20 \, a b c^{4}\right )} d + 12 \, {\left (15 \, b^{4} c^{2} - 100 \, a b^{2} c^{3} + 128 \, a^{2} c^{4}\right )} e - {\left (105 \, b^{5} c - 760 \, a b^{3} c^{2} + 1296 \, a^{2} b c^{3}\right )} f + 2 \, {\left (120 \, {\left (b^{2} c^{4} + 20 \, a c^{5}\right )} d - 12 \, {\left (5 \, b^{3} c^{3} - 28 \, a b c^{4}\right )} e + {\left (35 \, b^{4} c^{2} - 216 \, a b^{2} c^{3} + 240 \, a^{2} c^{4}\right )} f\right )} x\right )} \sqrt {c x^{2} + b x + a}}{15360 \, c^{5}}\right ] \]

input
integrate((c*x^2+b*x+a)^(3/2)*(f*x^2+e*x+d),x, algorithm="fricas")
 
output
[-1/30720*(15*(24*(b^4*c^2 - 8*a*b^2*c^3 + 16*a^2*c^4)*d - 12*(b^5*c - 8*a 
*b^3*c^2 + 16*a^2*b*c^3)*e + (7*b^6 - 60*a*b^4*c + 144*a^2*b^2*c^2 - 64*a^ 
3*c^3)*f)*sqrt(c)*log(-8*c^2*x^2 - 8*b*c*x - b^2 + 4*sqrt(c*x^2 + b*x + a) 
*(2*c*x + b)*sqrt(c) - 4*a*c) - 4*(1280*c^6*f*x^5 + 128*(12*c^6*e + 13*b*c 
^5*f)*x^4 + 16*(120*c^6*d + 132*b*c^5*e + (3*b^2*c^4 + 140*a*c^5)*f)*x^3 + 
 8*(360*b*c^5*d + 12*(b^2*c^4 + 32*a*c^5)*e - (7*b^3*c^3 - 36*a*b*c^4)*f)* 
x^2 - 120*(3*b^3*c^3 - 20*a*b*c^4)*d + 12*(15*b^4*c^2 - 100*a*b^2*c^3 + 12 
8*a^2*c^4)*e - (105*b^5*c - 760*a*b^3*c^2 + 1296*a^2*b*c^3)*f + 2*(120*(b^ 
2*c^4 + 20*a*c^5)*d - 12*(5*b^3*c^3 - 28*a*b*c^4)*e + (35*b^4*c^2 - 216*a* 
b^2*c^3 + 240*a^2*c^4)*f)*x)*sqrt(c*x^2 + b*x + a))/c^5, -1/15360*(15*(24* 
(b^4*c^2 - 8*a*b^2*c^3 + 16*a^2*c^4)*d - 12*(b^5*c - 8*a*b^3*c^2 + 16*a^2* 
b*c^3)*e + (7*b^6 - 60*a*b^4*c + 144*a^2*b^2*c^2 - 64*a^3*c^3)*f)*sqrt(-c) 
*arctan(1/2*sqrt(c*x^2 + b*x + a)*(2*c*x + b)*sqrt(-c)/(c^2*x^2 + b*c*x + 
a*c)) - 2*(1280*c^6*f*x^5 + 128*(12*c^6*e + 13*b*c^5*f)*x^4 + 16*(120*c^6* 
d + 132*b*c^5*e + (3*b^2*c^4 + 140*a*c^5)*f)*x^3 + 8*(360*b*c^5*d + 12*(b^ 
2*c^4 + 32*a*c^5)*e - (7*b^3*c^3 - 36*a*b*c^4)*f)*x^2 - 120*(3*b^3*c^3 - 2 
0*a*b*c^4)*d + 12*(15*b^4*c^2 - 100*a*b^2*c^3 + 128*a^2*c^4)*e - (105*b^5* 
c - 760*a*b^3*c^2 + 1296*a^2*b*c^3)*f + 2*(120*(b^2*c^4 + 20*a*c^5)*d - 12 
*(5*b^3*c^3 - 28*a*b*c^4)*e + (35*b^4*c^2 - 216*a*b^2*c^3 + 240*a^2*c^4)*f 
)*x)*sqrt(c*x^2 + b*x + a))/c^5]
 
3.2.5.6 Sympy [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 1360 vs. \(2 (230) = 460\).

Time = 0.63 (sec) , antiderivative size = 1360, normalized size of antiderivative = 5.76 \[ \int \left (a+b x+c x^2\right )^{3/2} \left (d+e x+f x^2\right ) \, dx=\text {Too large to display} \]

input
integrate((c*x**2+b*x+a)**(3/2)*(f*x**2+e*x+d),x)
 
output
Piecewise((sqrt(a + b*x + c*x**2)*(c*f*x**5/6 + x**4*(13*b*c*f/12 + c**2*e 
)/(5*c) + x**3*(7*a*c*f/6 + b**2*f + 2*b*c*e - 9*b*(13*b*c*f/12 + c**2*e)/ 
(10*c) + c**2*d)/(4*c) + x**2*(2*a*b*f + 2*a*c*e - 4*a*(13*b*c*f/12 + c**2 
*e)/(5*c) + b**2*e + 2*b*c*d - 7*b*(7*a*c*f/6 + b**2*f + 2*b*c*e - 9*b*(13 
*b*c*f/12 + c**2*e)/(10*c) + c**2*d)/(8*c))/(3*c) + x*(a**2*f + 2*a*b*e + 
2*a*c*d - 3*a*(7*a*c*f/6 + b**2*f + 2*b*c*e - 9*b*(13*b*c*f/12 + c**2*e)/( 
10*c) + c**2*d)/(4*c) + b**2*d - 5*b*(2*a*b*f + 2*a*c*e - 4*a*(13*b*c*f/12 
 + c**2*e)/(5*c) + b**2*e + 2*b*c*d - 7*b*(7*a*c*f/6 + b**2*f + 2*b*c*e - 
9*b*(13*b*c*f/12 + c**2*e)/(10*c) + c**2*d)/(8*c))/(6*c))/(2*c) + (a**2*e 
+ 2*a*b*d - 2*a*(2*a*b*f + 2*a*c*e - 4*a*(13*b*c*f/12 + c**2*e)/(5*c) + b* 
*2*e + 2*b*c*d - 7*b*(7*a*c*f/6 + b**2*f + 2*b*c*e - 9*b*(13*b*c*f/12 + c* 
*2*e)/(10*c) + c**2*d)/(8*c))/(3*c) - 3*b*(a**2*f + 2*a*b*e + 2*a*c*d - 3* 
a*(7*a*c*f/6 + b**2*f + 2*b*c*e - 9*b*(13*b*c*f/12 + c**2*e)/(10*c) + c**2 
*d)/(4*c) + b**2*d - 5*b*(2*a*b*f + 2*a*c*e - 4*a*(13*b*c*f/12 + c**2*e)/( 
5*c) + b**2*e + 2*b*c*d - 7*b*(7*a*c*f/6 + b**2*f + 2*b*c*e - 9*b*(13*b*c* 
f/12 + c**2*e)/(10*c) + c**2*d)/(8*c))/(6*c))/(4*c))/c) + (a**2*d - a*(a** 
2*f + 2*a*b*e + 2*a*c*d - 3*a*(7*a*c*f/6 + b**2*f + 2*b*c*e - 9*b*(13*b*c* 
f/12 + c**2*e)/(10*c) + c**2*d)/(4*c) + b**2*d - 5*b*(2*a*b*f + 2*a*c*e - 
4*a*(13*b*c*f/12 + c**2*e)/(5*c) + b**2*e + 2*b*c*d - 7*b*(7*a*c*f/6 + b** 
2*f + 2*b*c*e - 9*b*(13*b*c*f/12 + c**2*e)/(10*c) + c**2*d)/(8*c))/(6*c...
 
3.2.5.7 Maxima [F(-2)]

Exception generated. \[ \int \left (a+b x+c x^2\right )^{3/2} \left (d+e x+f x^2\right ) \, dx=\text {Exception raised: ValueError} \]

input
integrate((c*x^2+b*x+a)^(3/2)*(f*x^2+e*x+d),x, algorithm="maxima")
 
output
Exception raised: ValueError >> Computation failed since Maxima requested 
additional constraints; using the 'assume' command before evaluation *may* 
 help (example of legal syntax is 'assume(4*a*c-b^2>0)', see `assume?` for 
 more deta
 
3.2.5.8 Giac [A] (verification not implemented)

Time = 0.29 (sec) , antiderivative size = 403, normalized size of antiderivative = 1.71 \[ \int \left (a+b x+c x^2\right )^{3/2} \left (d+e x+f x^2\right ) \, dx=\frac {1}{7680} \, \sqrt {c x^{2} + b x + a} {\left (2 \, {\left (4 \, {\left (2 \, {\left (8 \, {\left (10 \, c f x + \frac {12 \, c^{6} e + 13 \, b c^{5} f}{c^{5}}\right )} x + \frac {120 \, c^{6} d + 132 \, b c^{5} e + 3 \, b^{2} c^{4} f + 140 \, a c^{5} f}{c^{5}}\right )} x + \frac {360 \, b c^{5} d + 12 \, b^{2} c^{4} e + 384 \, a c^{5} e - 7 \, b^{3} c^{3} f + 36 \, a b c^{4} f}{c^{5}}\right )} x + \frac {120 \, b^{2} c^{4} d + 2400 \, a c^{5} d - 60 \, b^{3} c^{3} e + 336 \, a b c^{4} e + 35 \, b^{4} c^{2} f - 216 \, a b^{2} c^{3} f + 240 \, a^{2} c^{4} f}{c^{5}}\right )} x - \frac {360 \, b^{3} c^{3} d - 2400 \, a b c^{4} d - 180 \, b^{4} c^{2} e + 1200 \, a b^{2} c^{3} e - 1536 \, a^{2} c^{4} e + 105 \, b^{5} c f - 760 \, a b^{3} c^{2} f + 1296 \, a^{2} b c^{3} f}{c^{5}}\right )} - \frac {{\left (24 \, b^{4} c^{2} d - 192 \, a b^{2} c^{3} d + 384 \, a^{2} c^{4} d - 12 \, b^{5} c e + 96 \, a b^{3} c^{2} e - 192 \, a^{2} b c^{3} e + 7 \, b^{6} f - 60 \, a b^{4} c f + 144 \, a^{2} b^{2} c^{2} f - 64 \, a^{3} c^{3} f\right )} \log \left ({\left | 2 \, {\left (\sqrt {c} x - \sqrt {c x^{2} + b x + a}\right )} \sqrt {c} + b \right |}\right )}{1024 \, c^{\frac {9}{2}}} \]

input
integrate((c*x^2+b*x+a)^(3/2)*(f*x^2+e*x+d),x, algorithm="giac")
 
output
1/7680*sqrt(c*x^2 + b*x + a)*(2*(4*(2*(8*(10*c*f*x + (12*c^6*e + 13*b*c^5* 
f)/c^5)*x + (120*c^6*d + 132*b*c^5*e + 3*b^2*c^4*f + 140*a*c^5*f)/c^5)*x + 
 (360*b*c^5*d + 12*b^2*c^4*e + 384*a*c^5*e - 7*b^3*c^3*f + 36*a*b*c^4*f)/c 
^5)*x + (120*b^2*c^4*d + 2400*a*c^5*d - 60*b^3*c^3*e + 336*a*b*c^4*e + 35* 
b^4*c^2*f - 216*a*b^2*c^3*f + 240*a^2*c^4*f)/c^5)*x - (360*b^3*c^3*d - 240 
0*a*b*c^4*d - 180*b^4*c^2*e + 1200*a*b^2*c^3*e - 1536*a^2*c^4*e + 105*b^5* 
c*f - 760*a*b^3*c^2*f + 1296*a^2*b*c^3*f)/c^5) - 1/1024*(24*b^4*c^2*d - 19 
2*a*b^2*c^3*d + 384*a^2*c^4*d - 12*b^5*c*e + 96*a*b^3*c^2*e - 192*a^2*b*c^ 
3*e + 7*b^6*f - 60*a*b^4*c*f + 144*a^2*b^2*c^2*f - 64*a^3*c^3*f)*log(abs(2 
*(sqrt(c)*x - sqrt(c*x^2 + b*x + a))*sqrt(c) + b))/c^(9/2)
 
3.2.5.9 Mupad [F(-1)]

Timed out. \[ \int \left (a+b x+c x^2\right )^{3/2} \left (d+e x+f x^2\right ) \, dx=\int {\left (c\,x^2+b\,x+a\right )}^{3/2}\,\left (f\,x^2+e\,x+d\right ) \,d x \]

input
int((a + b*x + c*x^2)^(3/2)*(d + e*x + f*x^2),x)
 
output
int((a + b*x + c*x^2)^(3/2)*(d + e*x + f*x^2), x)